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Abstract 

Applicability and effectiveness of robust control 
techniques to a single-stage-to-orbit (SSTO) airbreathing 
hypersonic vehicle on an ascent accelerating path and their 
effectiveness are explored in this paper. An SSTO control 
system design problem, requiring high accuracy tracking of 
velocity and altitude commands while limiting angle of 
attack oscillations, minimizing control power usage and 
stabilizing the vehicle all in the presence of atmospheric 
turbulence and uncertainty in the system, was formulated to 
compare results of the control designs using Hoo and p.- 
synthesis procedures. The math model, an integrated 
flight/propulsion dynamic model of a conical accelerator 
class vehicle, was linearized as the vehicle accelerated 
through Mach 8. Controller analysis was conducted using 
the singular value technique and the p-analysis approach. 
Analysis results were obtained in both the frequency and the 
time domains. The results clearly demonstrate the inherent 
advantages of the structured singular value framework for 
this class of problems. Since payload performance margins 
are so critical for the SSTO mission, it is crucial that 
adequate stability margins be provided without sacrificing 
any payload mass. 

1,0 Introduction 


The single-stage-to-orbit (SSTO) airbreathing 
hypersonic vehicles under consideration present significant 
challenges in many technological areas and especially in the 
realm of flight control. They traverse a broader flight 
envelope than any aircraft flown to date and must emphasize 
performance during the entire flight regime to achieve their 
mission objective. Large variations in vehicle static and 
dynamic characteristics and mass properties, such as 
significant movement of aerodynamic center of pressure, 
result in continuously changing static stability margins 
throughout the flight envelope Furthermore in general, 
an additional source of uncertainty arises from the accuracy 
of mathematical dynamic models used to describe the vehicle 
in control system design. These challenges and the limited 
availability of empirical data above Mach 8 in aerodynamics, 
propulsion, aeroelasticity, heating and on their combined 
effects on the vehicle's mission performance dictate the need 
for a robust yet performance oriented control system. 

The airframe/propulsion interactions, possibly the most 
complex of any vehicle, are of critical importance to the 
hypersonic vehicle mission success. The high sensitivity of 
the airbreathing propulsion system performance to the 
changes in angle of attack and dynamic pressure have been 
identified by Walton 2 and confirmed by Shaughnessy, et al 1 . 
Furthermore, atmospheric turbulence and especially large 
density variations at high altitude and Mach number 3 
introduce another significant source of uncertainty in the 
airbreathing propulsion system performance with which the 


control system must contend. 

In addition, as pointed out by Cribbs 4 , uncertainty in 
parameters, such as propulsive efficiency, drag and vehicle 
weight, all have major effect on vehicle performance 
margins in reaching orbital speed. The significantly 
detrimental effect of control surface deflection induced drag 


on the amount of fuel 5 to orbit provides another 
compelling reason for the importance of the control system 
optimization in hypersonic class vehicles. 

These issues and their impact on the control system 
development have been previously recognized by a number 
of researchers, among them Shaughnessy, et al. *, McRuer, 
et al. Anderson, et al. 2 and others. The control work in 
this area has primarily addressed the issue of an airbreathing 
hypersonic cruiser 2 which assumes equilibrium steady state 
flight with changes in coefficients of the equations of 
motion stemming from poor model description rather than 
changing flight parameters due to accelerated flight. 
Furthermore, the control laws developed for the the ascent 
phase have either disregarded the impact of angle of attack 
variations on airbreathing propulsion performance 2 or, 
while addressing tracking and atmospheric turbulence issues, 
did not explicitly consider performance robustness 8 . 

Recent application of modem robust control theory to 
the Space Shuttle 9 and fighter aircraft 10 flight control 
systems demonstrated potential benefits in dealing with the 
challenges mentioned above. The objective of this research 
is to assess the applicability and to exploit the capability of 
modem multivariable robust control theory to explicitly deal 
with both performance and uncertainty arising from 
changing flight conditions and vehicle characteristics. The 
problem is formulated to deal with the challenges associated 
with an SSTO hypersonic vehicle and its airbreathing 
propulsion system. A structured uncertainty model, 
representing parametric variations as actuator uncertainty, is 
used to compare two modem design procedures, Hoo and p- 
syn thesis. 

The paper consists of several sections. Section two 
provides a symbol list and section three gives a brief 
background review of robustness measures using 
conventional singular value methods. The issue of possible 
conservative solutions to some practical problems is 
discussed and the structured singular value (SSV) is 
introduced. The theoretical basis for the two synthesis 
procedures is also provided. This is followed in section four 
by a five state longitudinal linear model description and its 
derivation from an ascent trajectory of a conical accelerator 
class vehicle. Section five provides the problem description 
which discusses the uncertainty model, Hoo weighting 
function selection, and the explicit inclusion of stochastic 
atmospheric turbulence in a controller design. The last two 
sections deal with comparison of Hoo and p controller 
synthesis and analysis techniques as well as the conclusions 
derived from this study. 


2.0 Symbols 


SSTO 

single stage to orbit 

Hoc 

H infinity norm 

P 

structured singular value 

A 

uncertainty matrix 

w 

uncertainty matrix input vector 

e 

error/performance vector 

z 

uncertainty matrix output vector 

d 

exogenous inputs, i.e. noise, commands, 
turbulence vector 
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P(S) 

generalized plant structure 

K(s) 

controller 

G(s) 

augmented system plant, containing K(s) 
andP(s) 

F U (G, A) 

perturbed closed loop response 
between e and d 

e 1 

generalized output performance vector, 
e' = [we] 

d’ 

generalized exogenous input vector, 
d f = [ z d ] 

Fl(P, K) 

closed loop response between e' and d* 

u 

control effectors vector 

y 

sensed variables vector 

INI- 

infinity norm 

sup 

CO 

maximum value over frequency co 

a 

maximum singular value 

iff 

if and only if 

II-IIm 

structured singular value over all frequency 

CO 

degree of freedom 

DOF 

D 

positive definite Hermitian matrix 

V 

velocity, ft/sec 

a 

inertial angle of attack, deg 

q 

pitch rate, deg/sec 

0 

pitch angle, deg 

h 

altitude, ft 

8e 

symmetric elevon 

m/ 

fuel mass flow rate 

LTI 

linear time invariant 

W A 

uncertainty weighting matrix 

Wnoise 

measurement noise matrix 

Wp[«] 

performance weighting matrix for a given 
variable 

a 

bandwidth of elevon dynamics, 25 rad/sec 

b 

bandwidth of fuel flow rate dynamics, 
100 rad/sec 

Wj 

unity magnitude white Gaussian noise 

Fu,w 

longitudinal and vertical Dryden filters 

[V,h] c 

commanded velocity and altitude 

[V, h] ep 

performance weighted velocity and altitude 
error 

Hp 

performance weighted variables 

ASe, A8rh/ 

uncertainty in control effectors 

Se c , 8ifi/c 

controller commanded control inputs 

8e e ft, 8m/e« 

effective control inputs 

8e c , 8rh/ c 

uncertainty matrix inputs in the general 
structure 

X noise 

state measurement noise 

At 

time increment 


3.0 Theoretical Review of ji 

This section provides essential theorems for robustness 
and performance analysis in a control system with 
uncertainty. Far more detailed and rigorous discussion is 
presented in references [1 1-14]. 


Analysis methods based on singular values have been 
successful in providing multiloop extensions for classical 
single loop techniques * However, these methods are 
limited to providing exact results, i.e. necessary and 
sufficient conditions, for robust stability for systems with 
unstructured uncertainty, defined as norm-bounded but 
otherwise unknown perturbations. Consider, for example, 
the standard problem of analyzing a feedback system with 
simultaneously occurring multiplicative uncertainty at the 
plant input and output. In order to apply the singular value 
techniques, both perturbations must be reflected to a single 
location in the feedback loop, thus immediately inducing 
conservatism. However, since the combination of linear 
transformations is linear, any uncertainty occurring at 
several different locations in the feedback loop can be 
rearranged as a single block diagonal perturbation in a larger 
feedback loop. In other words, even unstructured uncertainty 
at the loop component level becomes highly structured at 
the system level *2 

The general framework for the problem is introduced in 
figure la. Any linear combination of inputs, outputs, 
commands, perturbations and controller can be arranged into 
the form in the diagram 



Figure la. General Structure 

Furthermore, the exogenous inputs d, the perturbation A and 
the output error e are normalized to 1 with all weightings 
and scalings absorbed into the generalized plant structure P. 
This arrangement results in unit invariant conditions for 
robust stability and performance expressed in terms of p. and 
presented later in this section. For the purposes of analysis 
the controller K can be considered an element of a larger 
plant G and be absorbed along with generalized plant P into 
its structure. The diagram for the analysis reduces to that in 
figure lb. 



Figure lb. Analysis Structure 
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The analysis problem itself involves determining whether 
the error e remains in a desired set for sets of inputs d and 
uncertainty A. The resulting structure of G can be 
partitioned as 


[wl 

Gii 

G12 

z 

L e J 

. G21 

G22 . 

.d. 


Closing the upper loop of G with the uncertainty matrix A 
results in a linear fractional transformation given by 

e = Fu(G, A) d = [ G22+ G21 A (I - Gi 1 A) 1 G12 ] d (3.2) 
The structure for synthesis is similarly given in figure lc. 



Figure lc. Synthesis Structure 

The inputs and outputs associated with A are absorbed into 
the exogenous input d' and error vector e'. The equation 
relating the generalized inputs to outputs is given by 

e = Fi (P, K) d = [ P11 + P12 K (I - P22 K)' 1 P21 ] d' (3.3) 

The conditions for stability and performance expressed 
in terms of Hoo bounds on portions of the generalized plant 
are given and discussed in some detail in references [12-14]. 
In the absence of uncertainty. A, the nominal performance 
objectives are expressed in terms of 

|| G22 ||~ = sup a ( G22 (ja>) ) ^1 ( 3 . 4 ) 

0) 

which relates the response e, || e ||2 ^ 1 , to the set of 
exogenous inputs d, || d ||2 ^ 1 . In practice, the use of 
palin gs and weightings is necessary to represent and 
normalize the varying frequency and spacial content of input 
and output sets. 

Consider plant perturbations that can destabilize a 
nominally stable system. Robust stability is satisfied for 
unstructured uncertainty iff the following condition holds, 

|| G 11 ||~ <1 for all A o (A) < 1 (3.5) 

In general for practical problems, the uncertainty 
consists of parameter variations and multiple norm-bounded 
perturbations that result from unmodelled dynamics of the 
system. Parameter variations often arise from changing 
flight conditions and represent changes in the coefficients of 
the equations describing the physical system. Unfortunately 


the norm-bounded test is insufficient and inadequate in 
dealing with robust performance and realistic models of plant 
uncertainty involving structure. 

To handle bounded structured uncertainty, the structured 
singular value (SSV) concept and the function p. are used to 
develop necessary and sufficient conditions The function 
p is defined as 

p (M) = — for M e C n * n 

mir^{o(A)|det(l + MA) = 0} ( 35 ) 

unless no A e A makes I + M A singular, in which case 
p ( M ) = 0 Based on the definition of p and its properties 
the robust stability condition is derived. Robust stability is 
satisfied iff 

IIG 11 ||n = supp [G 11 (j°>)] - 1 for all A e BA (3.6) 
to 

with prefix B denoting unit ball. Note that, in contrast with 
O, the value of p is dependent on G j 1 as well as on the 

structure of perturbations A. 

The question of interest to the control designer is how 
well does the system perform in the presence of uncertainty. 
Consider the issue of robust performance which describes 
performance with noise and perturbations occurring 
simultaneously: 

F u (G, A ) is stable and II F U (G, A) ||~<1 
for all A g BA 

iff ||G|| m <1 (3.7) 

The extreme cases for the structure of A provide the 
basis for the computational bounds on (X. The structure of 
A in general consists of repeated scalar blocks and full 
matrix blocks. |X can be computed exactly from the upper 
bound if the structure of A corresponds to 2S+F < 3, where 
S is the number of repeated scalar blocks and F is the 
number of full blocks, thus giving 

p(G) = D in1 D o( DGD ' 1 ) (3.8) 

In this particular problem A consisted of a single repeated 
scalar block given in a later section. The transformation 
DGD ' 1 is essentially a scaling of the inputs and outputs of 
G which does not change the value of |X. In addition, since 
[i can be computed exactly as a O plus scaling the methods 
developed for Hoo optimal control can be used to optimize \i . 

|X-analysis can be combined with Hoo optimal control to 
produce ^-synthesis which provides Hoo performance in the 
presence of structured uncertainty. The scaling matrices D 
and D ' 1 are used to reflect the structure of A over the 
frequency range. The problem now becomes reformulated as 
an Hoo-norm minimization of 

|| D Fi ( P, K ) D 1 ||~<1 ( 3 . 9 ) 
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known as D-K iteration. As the D-K name implies, the p- 
synthesis approach is to iterate between D and K until the 
solution converges The method is not guaranteed to 
produce a global minimum or to converge; however, the 
results, widely published in literature, have been successful 
in practical applications. 

4.0 Hypersonic Vehicle Model 

A conical accelerator configuration was used as an 
example for a generic airbreathing hypersonic vehicle 
As the vehicle accelerated through Mach 8 at 86,000 feet, a 
ten-state linear model representing the vehicle dynamics was 
obtained at this non-equilibrium flight condition, which is 
characterized by non-zero translational and rotational 
accelerations. The linear model was decoupled into a five- 
state longitudinal and a five-state lateral-directional model. 
The five state longitudinal vector x and control vector u, 
utilized in this study, are given by 


x =[ V a 

q 0 h ] T 

(4.1) 

u *{ 5 e 

5m/ ] T 

(4.2) 


The state and control variables are perturbation 
quantities and represent deviations from the nominal flight 
conditions. The open loop characteristics of the plant are 
unstable. Note that the altitude is included as a state 
variable to account for temperature, density and gravity 
gradients. These variations significantly affect the 
longitudinal long-period dynamics of the vehicle and add an 
aperiodic altitude mode caused by the variation of 
atmospheric density with altitude 

Several interesting nuances about this model are worth 
considering. In all practical problems, a linear time 
invariant (LTI) system is only an approximation to the real 
behavior. In most cases, time invariance of the system's 
characteristics is valid for "small" increments of time as 
determined by the researcher and dictated by the problem. In 
this particular case, the vehicle accelerates through Mach 8 
and the LTI system is valid for only that instant of time. 
However, if die parameter variation with time is represented 

as a multiplicative uncertainty such that G = G ( I + A ), then 
the linear system can be considered time invariant and can be 
used in LTI control design. 

Secondly, it is important to recognize that air mass 
flow rate changes almost instantaneously at hypersonic 
speeds with a change in angle of attack. This relationship is 
not explicitly included in the original nonlinear model; 
however, the control problem is formulated to account for 
the major airframe/propulsion interactions. This is reflected 
in explicit performance requirements on angle of attack as 

well as the engine control effector, rh/ . For generality and 

convenience, m/ is expressed in terms of fractional change 
from the nominal rather than in physical units. 

5.0 Problem Description 


The controller requirements were established based on 
the near optimal ascent trajectory and the sensitivity of the 
airbreathing propulsion system to angle of attack variations. 
Thus the specifications included high accuracy tracking of 
velocity and altitude commands while limiting angle of 
attack total deflection from nominal to 0.5 degrees, 
minimizing control power use and stabilizing the vehicle all 
in the presence of atmospheric turbulence and uncertainty in 
the system. It was assumed for the purposes of this study 
that Xstate = ymeas and was available for output feedback. 

The assumption was based on the availability of a inertial 
from calculations using outputs from an inertial 
measurement unit (IMU). 

The block diagram problem formulation is illustrated in 
figure 2a. All feedback state measurements were assumed to 
be corrupted by noise with the noise matrix represented by 

Wnoise = 1 0" 6 1 5 (5.1) 

The noise matrix was not intended to represent realistic 
sensor data and was included because the application of H°° 
to output feedback requires that the measured signals are 
corrupted by noise. Furthermore, the control actuator 
dynamics were represented by first order filters with 25 
rad/sec bandwith for elevon and 100 rad/sec bandwith for fuel 
flow rate. In a generic airbreathing hypersonic vehicle, the 
uncertainty, as discussed above, comes from different sources 
and occurs simultaneously. Thus, the very physics of the 
problem impose a structure on an uncertainty of a 
hypersonic vehicle. In order to explore the effects of 
structure on Hoo and p based analysis and synthesis 
techniques, actuator uncertainty was chosen for this 
problem. To represent parameter variation due to 
acceleration of the vehicle, a 20 percent multiplicative 
uncertainty. A, in the control effectiveness was introduced to 
the system where 


A = 


8 0 
0 8 . 


8eC 


(5.2) 


and the uncertainty weighting matrix was 


W A = 


0.2 

0 


0 

0.2 . 


(5.3) 


Thus, the control effectiveness was forced to vary from 80 
percent to 120 percent of the assumed nominal. At this 
stage of problem development, atmospheric density 
perturbations were assumed to be reflected in the 20 percent 
uncertainty in fuel flow rate effectiveness. 

Performance specifications for a flight control system 
translate quite well into an Hoo context for this problem. In 
designing for tight performance margins, performance 
weighting, or penalty functions, were augmented to the 
system. Hoo performance specifications make practical sense 
only when meaningful variables are specified for weighting 
functions. As illustrated in figure 2a, the weighting 
functions for control effector positions and rates as well as 
for the state vector were employed as output performance 
variables. 
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Specific performance requirements were derived from 
near fuel optimum ascent trajectory for the conical 
accelerator. The frequency dependent l^ormance 
weightings are used for velocity and altitude error. The time 

response requirements for velocity are specified ^10 ^rcent 

overshoot, translating into 0.5 high frequency gam, 40 sec 
time constant, giving 0.025 rad/sec as CTOSS^^ncy, 
and 5 percent steady state error, equivalent to low frequency 
gain of 50, all of which translate into 


... 0.5 ( s + 4.33e-2 ) 

WpV * * ~ ( s + 4.33e-4) 


(5.4) 


Similarly, altitude time response specifications 10 percent 
overshoot, 35 sec time constant, 5 percent steady state error, 

result in 


0.5(s + 4.95e-2J (5 .5) 

ph# ( s + 4.95e-4 ) 

Constant weighting was applied to a, q and 0 resulting in 



a 


15 ' 

Wp 

q 

= 

1 


. 0 . 


_ 1 - 


The weighting for a was based on the desire to attenuate 


atmospheric disturbances and to limit the magnitude of the 
outnuTas much as possible without violating performance 
Sments on other variables and was derived by itemtion. 
Similarly actuator position and rate limits were impo y 



(5.7) 

(5.8) 


Lastly, atmospheric turbulence is also explicitly 
included in the problem formulation. The primary concern 
is the effect of turbulence in engine performance. 

Turbulence can either affect performance directly by 
changing inlet flow conditions or by exciting actuator 
controllers leading to undesirable vehicle angle of attack 
variations Longitudinal and vertical Dryden turbulence 
filters are implemented as input weighting functions in the 
generalized framework (F u and F w in figure 2a) 

Filter inputs, wi and w 2 , to F u and F w , respectively 
represent white, zero mean, unit variance noise signals. The 
longitudinal Dryden filter output is assumed to act along 
velocity vector due to a small angle of attack flight 
SSL The vertical Dryden flier is S 

velocity, V 0 , to give angle of attack instead of vertical 
velocity perturbation. Furthermore, for correct 
implementation in a discrete simulation environment the 














continuous time Dryden filters, b y definition , must be 
divided by VNyquist frequency . 

• . P 1 ® b,ock diagram from figure 2a can be manipulated 
2h 0 *® g ?"^ f ^ l f w orico f figure la as depicted in figure 
b. Recall that all the input and output signals of the 
generalized plant P belong to the unity bounded sets with 
scaling absorbed into P. In this problem the performance 
weighting functions also served as the scaling factors for the 
output signal set. Thus, the input labels in figure 2b refer 
to the physical quantities represented by the inputs and 
outputs in figure 2a Note that the controller commanded 
inputs are also, after weighting, the A matrix inputs, thus 
defimng the system uncertainty. For the purpose of control 
synthesis the commands, atmospheric turbulence and 
actuator uncertainty are combined into the form given in 

in 8 !^ f C |i W,th the res . u,ts of the controller design discussed 
in the following section. 


A 5a 
ASmf 


V c h c 

* no toe 
W 1 W 2 

5e c 

Smfc 


P(S) 


K(s) 


5 -C 

#m»« 


^ep h ep 

“p Vp 

8 P 


Ve he 

a q e 



Figure 3a Hoo controUer frequency response analysis 

atmospheric turbulence in the longitudinal and vertical 
directions as can be seen in figure 3a The reader will recall 
_ om 860,1011 3 that nominal performance is satisfied iff 

° (F 22 ^) ) * 1 for all frequency. Satisfying the nominal 
performance condition signifies that specified response 

^SS nStlCS , ^ me i f0r the worse P° ssible combination of 
bounded inputs into the system with no uncertainty. The 

lack of robust performance and the contradictory results in 
robust stabihty between figures 3a and 3b, calculated with 
two different methods, will be addressed later. 


5E? 2h : ? ySiCal ? 686 ription of unity bounded input and 
output sets of generalized plant P. 

6.0 Design Comparison 

As the initial step in a controller design, an Hoo 
controller for a model with no uncertainty is obtained The 
nominal aircraft model is derived from the system 
interconnection, shown in figure 2b, by either deleting the 
rows and columns of P corresponding to w and z or setting 
A - 0- Frequency domain closed loop system analyses for 
nominal performance, robust stability and robust 
performance are illustrated in figures 3a-b. The closed loon 
system using the Hoo controller satisfies nominal 
performance requirements for simultaneous inputs of two 
tracking commands, velocity and altitude, and in the 
presence of 



0.01 0.1 1 
Frequency, rad/sec 


Figure 3b. H°o controller robust stability analysis 
maximum singular value 


using 


The time response of the nominal aircraft model to 

of 100 Wsec veloci, y change and 
v2ri™i J Mude change while encountering longitudinal and 
vertical atmospheric turbulence is presented in figures 4a-e. 

ese figures also include p controller nominal performance 
which wifi be discussed later in this section. Both velocity 

drived Tg 4a ' b ’ meet "* Performance requirements 
derived from tracking a near fuel optimum trajectory The 

effective angle of attack experienced by the propulsion 
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system, shown in figure 4c, also fulfills performance 
specifications of less than 0.5 degree total deviation. The 
elevon deflection, figure 4d, is around - 1 degree, thus 
limiting actuator induced drag. The fuel flow rate, figure 4e, 
avoids large sudden magnitude changes thus minimizing 

transients in the combustor. 

Minimizing the magnitude of an actuator deflection as 
well as the magnitude and rate of fuel flow rate in the engine 
improves vehicle performance. Typically, for an 
airbreathing hypersonic SSTO vehicle, payload fraction is 
only 3 percent while fuel fraction is around 60 percent . 
Therefore, any improvement in fuel fraction due to reduction 
in actuator induced drag has a potential to substantially 
increase payload fraction. Furthermore, since the 
performance of the propulsion system is extremely sensitive 
to changing conditions in the inlet, combustor and nozzle, it 
is important to minimize the transients in all parts of the 
engine. The control system minimizes perturbations in the 
inlet conditions by limiting angle of attack and aids 
combustion stability with smooth changes in fuel flow rate. 

As previously mentioned, the goal of the design is to 
maintain performance with 20 percent control power 
uncertainty present in the system. Prior to dealing with 
robust performance, the issue of robust stability must be 
addressed. Two methods for stability evaluation were 
applied to the aircraft model depicted in figure 2a. The 
results of the singular value analysis, a ( GnGco) ) ^ 1 > 
illustrated in figure 3b, indicate that the ctosed loop system 
is unstable for some A s with magnitude a (A s ) >1/16. 
However, the p-analysis technique, p(Gn(jco) ) < 1 , which 
explicitly considers the block diagonal structure of the 
uncertainty A, clearly demonstrates the stability of the same 
closed system in figure 3a. In fact, it becomes evident from 
the robust performance test that the system maintains 
stability and nominal performance for A with a much larger 
magnitude than the A s that would destabilize the system 
according to the singular value robustness test. 

Furthermore, time domain plots for a perturbed closed loop 
system with Hoc controller, figures 5a-g, confirm the 
stability of the system. The A used in the time simulations 
was the worst-case, real-rational, stable perturbation with 
a (A) > 1/1 -1 • The close loop poles of the perturbed 
system contain a complex pair that is essentially neutrally 
stable, -2.39e-1 0±3.91 , and manifests itself as a harmonic 

like oscillation on the time plots. 

Since the closed loop system is stable in the presence of 
structured perturbations defined in equation 5.2, the question 
of robust performance can now be addressed. It is evident 
from figure 3a that at least some performance requirements 
are no longer satisfied as indicated by the violation of the 
robust performance condition, i.e. p ( G(jto) ) > 1 • The 
effect of the uncertainty on velocity and altitude, illustrated 
in figures 5a-b, is negligible. The most apparent difference 
in time response caused by the perturbations is the elevon 
time history. The amplitude of the perturbed response is 

more than twice the nominal response, figure 5e. This 

increased elevon activity leads directly to the poor angle of 
attack response depicted in figure 5c. The fuel flow rate 
response is also degraded significantly as indicated by the 
oscillations seen in figure 5g. The overall degradation in 


performance is primarily a result of higher sensitivity to 
atmospheric turbulence. Since any elevon deflection induces 
significant drag, and fuel flow rate and angle of attack 
deviations cause a loss of propulsive performance, it would 
be highly beneficial to improve robust performance as much 
as possible in the presence of the given uncertainty. 

An H°° controller designed with 20 percent uncertainty 
explicitly included in the system does not satisfy nominal 
performance. Figure 6 provides a sample time response ot 
velocity to a 100 ft/sec step command. The response barely 
approaches 60 ft/sec velocity change which does not fulliU 
requirements of either steady state error or rise time. Further 
analysis indicates that either system uncertainty conditions 
or performance specifications on the tracking variables must 
be relaxed. Hence, an FU controller cannot fulfill robust 
performance requirements as specified for this problem. 

In an attempt to improve robust performance with 
original specifications, a p based controller is computed 
based on D-K iterations. The nominal performance of the p 
controller compares well with the nominal performance of 
Hoo controller as illustrated in figures 4a-e. However, it is 
in the robust performance that the advantage of a p 


controller becomes apparent. 

It is instructive to compare time response plots ot the 
perturbed closed loop systems with Hoo and p controllers as 
presented in figures 5a-g. The perturbed response of the p 
controller closed loop system is nearly the same as its 
nominal response. The performance improvements over an 
Hoo controller are precisely in the variables that the former 
had trouble handling. The angle of attack demonstrates 
almost negligible response to perturbations as illustrated in 
figure 5d. The elevon response remains essentially 
unaffected by 20 percent uncertainty, figure 5f, which is also 
true for fuel flow rate, figure 5g. Furthermore, reducing fuel 
flow rate bandwidth from 100 rad/sec to 10 rad/sec to 
account for potential delays upstream of the controlling 
valves has negligible effect on the control system response. 

The frequency domain analysis supports the conclusion 
of satisfied performance requirements in the presence of 
control effector uncertainty. The robust performance 
condition in figure 7 is met. The robust stability test 
results are identical for both singular value and p techniques. 
This is due to the fact that the p controller reflects die 
uncertainty structure and, thus, it is absorbed into the close 
loop system. The p based controller successfuUy handles 
actuator uncertainty without sacrificing system performance. 
The encouraging preliminary results of the p controller 
establish the technique as potentially successful in dealing 
with unique characteristics of hypersonic class vehicles. 
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Figure 4a. Velocity response to velocity command, altitude 
command, vertical and longitudinal turbulence. 


Figure 4d. Elevon response to velocity command, altitude 
command, vertical and longitudinal turbulence. 




Figure 4b. Altitude response to velocity command, altitude c - 

command, vertical and longitudinal turbulence. Fl 8 ure 4e - Fuel flow rat e response to velocity command 

altitude command, vertical and longitudinal turbulence 




Figure 4c. Angle of attack response to velocity command c - . 

altitude command, vertical and longitudinal tuibulence. gUre 5 * Ve,oclt y response of perturbed system to velocity 

command, altitude command, vertical and longitudinal 
turbulence. 
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time, sec 


Figure 5b. Altitude response of perturbed system to velocity Figure 5e. Elevon response of perturbed system to velocity 
command, altitude command, vertical and longitudinal command, altitude command, vertical and longitudinal 

turbulence. turbulence. 



Tim#, mc 


Figure 5c. Angle of attack response of perturbed system to 
velocity command, altitude command, vertical and longitudinal 
turbulence. 



Tim*, mc 

Figure 5f. Elevon response of perturbed system to velocity 
command, altitude command, vertical and longitudinal 
turbulence. 



Figure 5d. Angle of attack response of perturbed system to Figure 5g. Fuel flow rate response of perturbed system to 
velocity command, altitude command, vertical and longitudinal velocity command, altitude command, vertical and longitudinal 
turbulence. turbulence. 
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Figure 6. Velocity response to velocity command, altitude 
command, vertical and longitudinal turbulence of a system 
with an H°° controller designed explicitly with uncertainty. 



Frequency, rad/sec 


Figure 7. controller frequency response analysis 


7.0 Conclusions 

Applicability of robust control techniques to an SSTO 
airbreathing hypersonic vehicle on an ascent accelerating 
path and their effectiveness are explored in this paper. 
Several important issues related to control system design 
should be noted. The most important results are based on 
comparison between H«, and p techniques. Since 
uncertainty plays an integral part of hypersonic vehicle 
characteristics, its effect on the analysis and synthesis of 
various control system design techniques is important to 
understand. 

The characteristics describing airbreathing hypersonic 
vehicles and the requirements imposed on the control system 
translate explicitly into H*. domain specifications as 
illustrated in this paper. However, the H«, controller suffers 
performance degradation with introduction of control effector 
uncertainty into the system. The (i controller preserves the 


required performance while providing stability robustness. 
Hence, p synthesis, by taking into consideration the 
structure of the uncertainty in this problem, results in an 
improved robust performance over the Hoc controller. As 
previously mentioned the extent of conservatism induced by 
singular value robustness analysis varies considerably from 
problem to problem. In this typical hypersonic problem 
formulation, singular value analysis is conservative for a 
relatively benign control effector uncertainty. 

The cumulative results of this research imply the 
importance of p as both the analysis and the synthesis tool 
for an airbreathing hypersonic vehicle. Since uncertainty 
occurs simultaneously from many different sources and the 
degree of uncertainty is high, the physical behavior of the 
system introduces structure into the problem. It is essential 
to employ a methodology that takes full advantage of these 
physical characteristics. The p-analysis and synthesis 
technique preserves the structural relationship between 
uncertainty and performance variables, allowing the designer 
a systematic approach to explore tradeoffs between the two. 
Failure to account for this structural relationship can result 
in excessively conservative specifications and poor designs 
for an airbreathing hypersonic vehicle. 
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